Group projection method in statistical systems 



Krzysztof Redlich " Frithjof Karsch and Ahmed Tounsi 

" Fakultat fur Physik, Universitat Bielefeld, Postfach 100 131, D-33501 Bielefeld, Germany 
^ Institute of Theoretical Physics, University of Wroclaw, PL-50204 Wroclaw, Poland 
Laboratoire de Physique Theorique et Hautes Energies, Universite Paris 7, Paris, France 

Abstract. We discuss an application of group theoretical methods to the formulation of the 
thermodynamics of systems constrained by the conservation laws described by a semi-simple 
compact Lie group. A general projection method that allows to construct a partition function 
for a given irreducible representation of the Lie group is outlined. Applications of the method 
in Lattice Gauge Theory (LGT) for non-zero baryon number and in the phenomenological 
description of particle production in ultrarelativistic heavy ion collisions are also indicated. 



1. Introduction 

In the formulation of the thermodynamics of a strongly interacting medium one needs in 
general to implement the constraints imposed by the conservation laws that are governed by 
an internal symmetry of the Hamiltonian |ITJIHI11I11I3|51|71[HI|. In this paper, we discuss how 
applying the basic properties of the Lie groups and their representations one can derive the 
partition function that accounts for these constraints. We present examples of the application 
of the method for different statistical systems and phenomenological models. 

2. Projected partition function 

The usual way of treating the problem of quantum number conservation in statistical physics 
is by introducing the grand canonical partition function, 



Z{^is,T,V) = Ti[e-^^^-^^% (1) 

with H being the Hamiltonian, S the charge operator, P the inverse temperature, V the volume 
of the system and f^s the chemical potential associated with the conserved charge S. The 
chemical potential in Eq. dU) plays the role of the Lagrange multiplier which is fixed by the 
condition that the charge is conserved on the average and has the required value (S) such that: 

dfis 

The grand canonical partition function dU) provides an adequate description of the 
statistical properties of the system only if the number of particles carrying charge S is 
asymptotically large and if its fluctuations can be neglected. 

In order to derive a more general statistical operator that is free from the above 
requirements one usually replaces the function ([T]) by the canonical partition function Zs 
that accounts for an exact charge conservation 

ZsiT,V):=Trs[e-^% (3) 
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The subscript S under the trace indicates that it is restricted to the states that carry an exact 
value S of the conserved charge. Obviously ZsiT, V) from Eq. © and Z(ns,T, V) from 
Eq. (HJ) are connected via cluster decomposition 



s=+oo 

Z{f^s,T,V)= J2 ZsXh (4) 

s=— oo 

with the fugacity parameter Xs := e^^^ . Thus, Zs is viewed as a coefficient in the Laurent 
series of Z in the fugacity. Applying the Cauchy formula in Eq. ^ one calculates Zs by 
taking an inverse transformation: 



Zs{T,V)= [ ^/e-^^*Z(0,T,l^). (5) 

The generating function Z{(j),T,V) := Z{\s = e^'^,T,V) is obtained from the grand 
canonical partition function ^ by "Wick rotation" of the chemical potential, fis i4>. The 
function Z is unique for all canonical partition functions that are formulated for a fixed value 
of the conserved charge. 

The integral ® describes the projection onto the canonical partition function that 
accounts for an exact conservation of an abelian charge. This is the projection procedure 
as Zs is also obtained from 

Zs{T,V) = TT[e-^^Ps] (6) 

where Ps is the projection operator on the states with a given value of the conserved charge S. 
For an additive quantum numbers, Ps is just a delta function Ps = 5s.,si- Using the Fourier 
expression of Ss.,s in Eq. © one reproduces the projected formula 

The conservation of an additive quantum number is usually related with an invariance 
of the Hamiltonian under the abelian U(l) internal symmetry group. In many physics 
applications it is of importance to generalize the projection method to symmetries that are 
related with a non-abelian Lie groups G. An example is a special unitary group SU(N) that 
plays an essential role in the strong interaction. A generalization of the projection method 
would require to specify the projection operator or the generating function. Consequently, the 
partition function obtained with a specific eigenvalue of the Casimir operator that fixes the 
multiplet of the irreducible representation of the symmetry group G can be determined. 

To find the generating function for the canonical partition function with respect to the 
symmetry group G, one introduces Z{g) \JJ 



~Z{g) := Tr[f/(^7)e-^^] (7) 

as the function on the group G with U(g) being a unitary representation of the group and 
g E G. The U(g) can be decomposed into irreducible representations Ua{g), 



U{g)=Y.U^{g), (8) 

a 

and thus, from Eqs. Q and ^ one can write explicitly 

a 

= E E {^-,^a I U^{g)e-''^ I (9) 
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where z/q, labels the states within the representation a and are degeneracy parameters 
of a given representation. Due to the requirement of an exact symmetry the only non- 
vanishing matrix elements of the evolution operator e~l^^ are those diagonal in u^- The 
matrix elements of Ua{g) are non-zero if they are diagonal in ^q. Finally, the matrix elements 
of the Hamiltonian are independent of the states within representation and those of U (g) of 
degeneracy factors. Consequently, the matrix elements in Eq. Q factorize and the generating 
function is 



Z{g)=j:^Z^iT,V), (10) 

where d{a) is the dimension and Xa is the character of the irreducible representation a. In 
the above expression Za is introduced as a canonical partition function with respect to G 
symmetry and is defined as 



Eq. (fTOb connects a canonical partition function Za to the generating function on the 
group. Thus, Za is the coefficient in the cluster decomposition of the generating function with 
respect to the characters of the representations associated with the symmetry group. 
The orthogonality relation of the characters, 

J dfl{g)Xa{g)X-/{g) = ^a,7 (12) 

allows to find the canonical partition function. From (fTH and (fT^ one gets 

Za{T, V) = d{a) J dMx*a{g)Z{g). (13) 

This result is a generalization of ^ to an arbitrary internal symmetry group that is a compact 
Lie group. The formula holds for any dynamical system as it is independent of the specific 
form of the Hamiltonian. 

To find the canonical partition function one needs to determine first the generating 
function Z(g) defined on the symmetry group G. If the group is of rank r, then the character 
of any irreducible representation is a function of r variables {71,.... ,7r}, thus also is the 
generating function (10). Diagonalizing the unitary operators U{g) under the trace in Eq.© 
and denoting by (A; = 1, r) the commuting generators of G, the generating function can 
be formulated on the maximal abelian subgroup of G as 

Z(7i,..,7,) = Tr[e-^^+*SU7.J.]. (14) 

Thus 2^(7) is just the GC partition function with complex chemical potentials that are 
associated with all generators of the Cartan sub-algebra. 

Equations (1131 and (fT4b provide a complete description of the statistical operators that 
account for the constraints imposed by the conservation laws of internal symmetries of the 
Hamiltonian. The simplicity of the projection formula (fT^l is that the operators that appear 
in the generating function are additive. One sees that the problem of extracting the canonical 
partition function with respect to an arbitrary compact Lie group G is reduced to the projection 
onto a maximal abelian subgroup of G. 

The generating function (fT4b can be calculated by applying standard perturbative 
diagrammatic methods or a mean field approach. However, if the interactions in the 
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Hamiltonian can be omitted or can effectively be described as a modification of the particle 
dispersion relations by implementing an effective particle mass, then the trace in Eq. ([T4b can 
be done exactly. In this particular situation the generating function can be written [IJ as 



^W-xpE^Z-] (.5) 
where the one-particle partition function in Boltzmann approximation 

= ^ / ^ exp {-P^f + mi) (16) 

is just a thermal phase-space that is available to all particles of mass that belong to a given 
irreducible multiplet. The sum is taken over all particle representations that are constituents 
of a thermodynamical system. 



2.1. Phenomenological model of colour confinement 

To illustrate how the projection method results in the partition function, we discuss a statistical 
model that accounts for the conservation of a non-abelian charge related with the global 
G = SUc(^)xUBil) internal symmetry of the Hamiltonian. As a physical system one 
considers a thermal fireball that is composed of quarks and gluons that carry the colour degrees 
of freedom related with SUc{N) symmetry [3, 4| and baryon number related with Usil) 
subgroup. The system has temperature T and volume V. The interactions between quarks 
and gluons are implemented effectively as resulting in dynamical particle masses that are T 
dependent, e.g. through (g) ~ gT. Thus, since under these conditions the free-particle 
dispersion relations are preserved, Eq. ( f]3l provides a correct description of the generating 
function on the symmetry group G. The sum in the exponent in dTSt gets the contributions 
from quarks, anti-quarks and gluons that transform respectively under the fundamental (0,1), 
its conjugate (1,0) and adjoint (1,1) representation of the SUc{N) symmetry group. Thus 0, 

* 

In Z{T, V, 7, 7b) = + ^4 + (17) 

"Q "Q «G 

where 7 = (71, ..,7Ar_i) are the parameters of the SUc{N) and 7^ of the Usil) symmetry 
groups. Through an explicit calculation of one-particle partition functions for massive quarks 
the corresponding generating functions read: 

lnZQ(r,V^,7,7.) = f^t^-^K,inm,/T) 

[ e^"-"/^XQ(^7) + e-*^-"/^XQ(n7)] (18) 

where the two terms in the brackets represent the contribution of quarks and antiquarks 
respectively. The corresponding contribution for massive gluons is obtained as 

lnZG(T,\/,7,7B) = E -K2{nmG/T)[xG{n^) + xU^l)] (19) 

The coefficients Qg^Qq and dq = N.dc = N"^ — 1, are respectively, the quark and gluon 
spin-isospin degeneracy factors and dimensions of the representations. 

The statistical operator is now obtained from (13) through the projection of the generating 
functions (fTSt and (fT9l) onto a given sector of quantum numbers. Here we quote the results 
for the SUc{'i) colour singlet partition function that represents a global colour neutrality 
condition (phenomenological confinement) of the quark-gluon plasma droplet [3|: 
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Figure 1. Thermodynami- 
cal functions obtained from the 
projected partition function (I20> 
are compared to their Stefan- 
Boltmann hmit as a function of 
VT^. The results for the energy 
density e/esB and the expectation 
value of the character in the fun- 
damental representation of SU (3) 
are indicated in the figure. The re- 
sults are obtained for a vanishing 
baryon-chemical potential /is in 
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Za,=oifiB,T,V) = / rf/x(7i,72)exp{ciXG + 



C2[LRCosh{PfiB) + iLjsmh{(3^B)]}- (20) 



where a finite average value of tlie baryon number is controlled by the chemical potential fiB- 
The constants Ci and C2 can be extracted from Eqs. (fT8l - fT9b and taking only the first terms 
in the series. The parameters Lj, Lr denotes the imaginary and real part of the character 
in the fundamental representation of the SUc{'i) group whereas xg is the character of the 
adjoint representation. The integration is done on the SUc{'i) group with an appropriate Haar 
measure /i(7). 

To illustrate how the group projection influence the thermodynamics of the system we 
display in Fig. 1 the behavior of the energy density and the thermal average of the characters 
as obtained from the colour singlet statistical operator from Eq. (l20b . The results are shown 
as a function of dimensionless parameter VT"^. In the absence of group projection and 
for a vanishing value of fiB the partition function would have a simple form; Z(T, V) = 
exp [ci + C2] . The resulting energy density b describes a well known Stefan-Boltzmann limit 
and (xg) = 1 in this case. For an asymptotically large value of VT^ — > 00, that corresponds to 
the thermodynamical limit, the colour projected results coincide with their Boltzmann values. 
However, for small values of VT^ a large suppression of thermodynamical quantities is seen 
in Fig. 1 . This is a generic feature of the partition function restricted to a fixed representation. 
The requirement of an exact conservation of the quantum numbers impose a strong constraint 
on the particle thermal phase-space that results in the suppression seen in Fig. 1. In the 
thermodynamical limit, the number of particles that carry a quantum number related with a 
given representation is so large that the above constraints are irrelevant. 

In finite temperature gauge theory the zero component of the gauge field Aq takes on 
the role of the Lagrange multiplier which guaranties that all states satisfy the Gauss law 
[TDJ [m. In the Euclidean space one can choose a gauge in such a way that Aq{x, t)\^ 
is a constant in Euclidean time, so that = g~^a"'^5ab- In such a gauge the Wilson loop 
defined as 



represents the character of the fundamental representation of the SUc{N) group Q. Thus, the 
projected partition function could be viewed as an effective model that connects the coloured 
quasi-particles L9J with the Wilson loop 1.11.1 . This partition function can also be related to 




(21) 
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Figure 2. Heavy quark potential 
from Ref. as a function 

of quark anti-quark separation 
RT. The results were obtained 
from Monte-Carlo simulations of 
quenched QCD at fixed value of 
the coupling constant (3 = 5.62 
and baryon number B = Q and 
B = 6. 
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the strong coupling effective free energy of the lattice gauge theory with a finite chemical 
potential L3J and to the effective potential of the SU (N) spin model for the Wilson loop. 

2.2. Projected partition function in Lattice Gauge Theory 

The formulation of a finite baryon number density in QCD on the lattice leads to a sever 
numerical problem of complex probability. The QCD partition function formulated on the 
Euclidean lattice and after integrating out the Wilson fermions, becomes 1 12>.13J 



where Sq is the gluon action, K := 1 — kM{U, /i^) is the fermion matrix and k is the hopping 
parameter. In order to perform the Monte-Carlo simulations with this statistical operator it 
is necessary that the measure be real and positive [.13t .1?! [T51l . One situation where det[K] 
is real is when there exist an invertible operator P such that K+ = PKP~^. For the Wilson 
fermions at zero chemical potential and for complex hb this relation hold with P = 75. 
However, for real and non- vanishing /is 7^ this is not anymore valid and the probabilistic 
interpretation of the path integral representation of the QCD partition function is violated 
|[13, 16 1 . Consequently the Monte-Carlo method is not valid and the numerical solution 
of the QCD thermodynamics is not accessible anymore. J One way to partly overcome the 
above problems in finite density QCD is to use a projection method [T/TOIIEI. Rather than 
introducing a non-vanishing chemical potential [14|, that is formulate QCD in the grand 
canonical ensemble with respect to U(l)-baryon symmetry, one my go over to a canonical 
formulation of the thermodynamics and fix directly the baryon number |l3l|T0l|T6|. Following 
the general discussion of Section 1 the projected partition function onto a given sector of 
baryon number could be obtained from 1.3. .161 : 



where the function Z under the integral represents the grand canonical partition function 
in Eq.(22) calculated with a complex fXB, and consequently is a real and positive quantity. 
Nevertheless, due to a Fourier integration the projected partition function still suffers a 
numerical problem of oscillating functions. However, in the quenched limit of QCD and 

I The complex structure of the fermionic contribution to the partition function is very transparent in the model 
discussed in the last section. From Eq. ilOi it is clear that the complex structure appears for iib ^ since 
L/ 7^ for the SU{N) system with N > 2. On the other hand the quark generating function (18) is indeed real 
for the complex baryon chemical potential. 




(22) 




(23) 
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Figure 3. The left hand figure: the statistical model results (broken line) on the energy 
dependence of A/tt ratio. The data are from |21 1. The right hand figure: a comparisons of the 
chemical freeze-out parameters with the phase boundary line. The upper thin line represents 
the LGT results obtained in Ref. |22| and the lower thin line describes the conditions of 
constant energy density that was fixed at = [22 1. The full ( Fodor-Katz) point is the end- 
point of the crossover transition from Ref. 1123 1 . The lower line is the unified freeze-out curve 
from L20J . 



for moderate values of baryon number the group integration was done analytically and the 
first results on QCD thermodynamics at finite baryon density were established fl6l. 

The presence of a finite net baryon number modifies non-trivially the properties of QCD 
medium. As an example of Monte-Carlo study we show in Fig. 2 a heavy quark anti-quark 
potential obtained from the Polyakov loop correlations [ 16l. For zero baryon number it shows 
the usual linearly rising behavior for the quenched case. For finite B = Q the potential 
remains finite at large distances due to screening of the static quark anti-quark pair by the 
already present net baryon charge in the system. 



2.3. Canonical projection in heavy ion collisions 

Central heavy ion collisions at relativistic incident energies represent an ideal tool to study 
nuclear matter at high temperatures and densities. Particle production is - at all incident 
energies - a key quantity to extract information on the properties of nuclear matter under 
these extreme conditions. In this context a particular role has been attributed to particles 
carrying strangeness that is related with U(l) invariance of the strong interactions ||2l|6|. The 
production of secondaries measured in heavy ion collisions was shown in the literature to be 
very satisfactory described in the context of statistical thermal models ll^lTl fTTlfTHlfT^Em . 
However, already the first attempts to describe strange particle production in low energy 
central and high energy peripheral heavy ion collisions have shown that the conservation of 
strangeness should be implemented exactly 

An exact formulation of strangeness conservation is implemented through the projection 
of the partition function onto states of fixed representation of the U(l) group lHJ- For the 
hadron resonance gas the strangeness neutral partition function is obtained from Q: 

Z'i^,{T, V) = ^ d<P exp ^n^*"'^ j (24) 

where S'„ is a thermal phase space available to all particles and resonances that carry 
strangeness n with n E (—3, 3). The density n], of particle i carrying strangeness s is derived 
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form (24) as Q 

{nlf = {n:f^F,{T,V) (25) 

where (n*)*^"" is the grand canonical density and Fs{T^V) is the suppression factor that 
measures a deviation of {n\)^ from its asymptotic, grand canonical value. For large V and/or 
T the factor Fs(T, V^) — > 1. The group projection implies that F < 1, thus it suppresses a 
particle densities that carries the U(l) charge. This suppression was found to increase with 
decreasing collision energy, increasing strangeness content of the particles and decreasing 
centrality of the collisions |7|. These properties are well observed in experimental data. As 
an illustration for the model comparison with experimental data [21j we show in Fig. (3-left) 
the results on the ratio of lambda to pion multiplicities |7 1. 

The statistical model that accounts for an exact conservation of U(l) charges was found 
in the literature to reproduce the basic features of particle yields obtained in heavy ion 
and hadron-hadron collisions. The yields were found to be well reproduced with thermal 
parameters (the temperature and baryon chemical potential) that follow a universal freeze-out 
line oi fixed energy per particle of 1 GeV [20|. The freeze-out line is shown in Fig. (3-right) 
together with the most recent LGT results on the position of the critical curve in the (T-hb) 
plane [22, 23|. 
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